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TIP  FLOWS  FOR  WINGS  AND  PROPELLERS 
AND  THEIR  EFFECT  ON 
THE  PREDICTED  PERFORMANCE 


Spyros  A.  Kinnasxand  Sangwoo  Pyo2 


ABSTRACT 

Various  wake  alignment  techniques  which  have  been  developed  in  the  past  for  the  analysis  of  non - 
cavitating  propellers  in  uniform  inflows  are  reviewed  and  some  of  them  are  extended  in  the  case 
of  cavitating  flows  subject  to  inclined  inflows .  The  effect  of  the  inclined  trailing  wake  geometry 
on  the  predicted  cavities  and  blade  forces  is  found  to  be  significant .  The  effect  of  the  tip  vortex 
detachment  location  on  the  shape  of  the  trailing  wake  and  on  the  pressure  distribution  on  the 
tip  is  studied  for  wings  and  propeller  blades.  The  local  viscous  flow  inside  of  the  core  of  a  tip 
vortex  is  formulated  via  a  parabolic  Navier-Stokes  approach.  Predicted  open  flow  characteristics 
and  unsteady  forces  acting  on  the  blades  of  an  inclined  shaft  propeller  are  compared  t#  those 
predicted  by  other  methods,  as  well  as  those  measured  in  experiments. 


NOMENCLATURE 

Vjn  ;  inflow  velocity  at  the  propeller  plane 

Cv  ;  pressure  coefficient  :  =  (p  —  poc)/(l/2)pV^ 

Assistant  Professor,  Ocean  Engineering  Group,  Department  of  Civil  Engineering,  The  University  of  Texas  at  Austin, 
kinnas@mailutexas.edu 

Postdoctoral  associate,  Ocean  Engineering  Group,  Department  of  Civil  Engineering,  The  University  of  Texas  at  Austin, 


D 

;  propeller  diameter 

F„ 

;  Froude  number  based  on  propeller  RPM;  = 

n2D/g 

Js 

;  advance  coefficient  based  on  ship  speed  :  = 

Vs/nD 

Kt 

;  trust  coefficient  :  =  T/pn2D 4 

Kq 

;  torque  coefficient  :  =  Q/pn2Ds 

n 

;  propeller  RPM 

P 

;  pressure 

Pv 

;  vapor  pressure 

Poo 

;  pressure  at  infinity  (at  the  shaft  depth) 

Vs 

;  ship  speed 

e 

;angular  coordinate  of  point  on  trailing  wake  surface 

01 s 

;blade  angle 

<?n 

;  cavitation  number  based  on  propeller  RPM 

;  =  (Poo  -pv)/(l/2)pn2D2 

U> 

;  propeller  angular  velocity 

Acronyms  -  Names  of  computer  codes 
FLAG  ; 

BEM  ; 

VLM  ; 

PSF-2  ; 

PSF-2IS  ; 

PUF-2  ; 

PUF-3  ; 

PUF-3A  ; 

HPUF-3A  ; 

HPUF-3AL/INC; 

PSF-10/PUF-10  ; 


FLow  Adapted  Grid 
Boundary  Element  Method 
Vortex  Lattice  Method 

Propeller  blade  analysis  code  in  Steady  Flow  based  on  VLM 
PSF-2,  1979  version  includes  Inclined  Shaft  effects 
Propeller  blade  analysis  code  in  Unsteady  Flow  based  on  VLM 

Propeller  U  nsteady  cavitating  Flow  analysis  code  based  on  VLM,  1979  original  version 
PUF-3,  1986-93  versions  includes  leading  edge  correction 
PUF-3,  1994  version  includes  effects  of  Hub 

PUF-3,  1996  version  includes  effects  of  wake  ALignment  and  effects  of  INClined  flow 
Propeller  blade  analysis  code  in  Steady/  Unsteady  Flow  based  on  BEM 
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1  INTRODUCTION 


The  hydrodynamic  analysis  of  lifting  surfaces  (wings  or  propeller  blades)  via  potential  flow  methods  requires 
knowledge  of  the  shape  of  the  trailing  wake  vorticity  sheet.  The  location  of  the  vortex  sheet  can  be  determined 
from  the  condition  of  zero  pressure  jump  across  the  sheet,  also  known  as  the  force- free  wake  condition.  This 
condition  is  equivalent  to  requiring  the  potential  jump  in  the  wake  to  be  convected  with  the  resulting  flow.  The 
prediction  of  the  evolution  of  a  wake  sheet  has  been  studied  extensively  in  the  past  and  a  review  of  the  related 
bibliography  with  emphasis  on  marine  propellers  may  be  found  in  (Pyo  1995)  and  (Pyo  &  Kinnas  1997). 

The  effect  of  the  main  parameters  of  the  wake  geometry,  such  as  pitch,  transition  wake  length,  and  ultimate 
wake  radius,  on  the  predicted  propeller  thrust  and  torque  was  first  studied  in  (Kerwin  <k  Lee  1978).  This  wake 
model  was  also  applied  when  the  effects  of  unsteady  sheet  cavitation  were  included  (Lee  1979).  The  importance 
of  the  wake  geometry  on  the  blade  forces  was  also  studied  by  (Hoshino  1991)  who  determined  the  wake  shape 
experimentally  by  using  an  LDV  system.  The  effect  of  the  wake  location  becomes  more  crucial  in  the  case  of  low 
advance  ratios,  as  for  example  in  the  case  of  helicopter  blades  (Morino  &  Bharadvaj  1988),  due  to  the  effect  of 
each  blade’s  wake  on  the  other  blades.  The  effect  of  the  wake  location  and  detailed  structure  is  also  important 
in  the  case  of  a  stage  of  a  multi-component  propulsor  working  in  the  wake  of  another. 

In  (Greeley  &;  Kerwin  1982)  an  iterative  procedure  was  developed  to  determine  the  wake  geometry  of  a 
propeller  in  axisymmetric  inflow.  Later  in  (Keenan  1989)  a  method  was  developed  for  the  treatment  of  the 
trailing  wake  in  the  case  of  non-axisymmetric  inflows.  These  inflows  are  mainly  due  to  the  ship’s  viscous  wake 
at  the  stern.  However,  the  propeller  blades  will  also  be  subject  to  unsteady  flows  in  the  case  of  inclined  inflows. 
This  inclination  may  be  either  due  to  an  inclined  shaft  or  due  to  a  transverse  inflow  during  maneuvering  of  the 
ocean  vehicle.  The  blades  will  be  subject  to  once  per  revolution  sinusoidally  varying  inflow  and  thus  unsteady 
forces,  which  may  lead  to  fatigue  and  eventually  deformation  or  damage  of  the  blade  (Boswell  et  al  1984).  The 
effect  of  the  inclined  wake  geometry  on  the  prediction  of  the  unsteady  blade  forces  in  the  case  of  inclined  inflow 
was  studied  in  (Kerwin  1979)  and  found  to  be  very  critical.  The  effect  of  the  inclined  wake  was  evaluated  at  0, 
90,  180,  and  270  degrees  and  interpolated  in  between. 

In  the  present  paper  first  a  review  of  the  existing  wake  models  for  marine  propeller  applications  is  given. 
The  model  of  (Greeley  &  Kerwin  1982)  is  then  extended  in  the  case  of  cavitating  flow.  Finally  an  inclined  wake 
model  is  developed  in  which  the  trailing  wake  is  aligned  (at  each  blade  angle)  with  the  axisymmetric  flow  to 
which  the  transverse  flow  due  to  the  shaft  inclination  is  superimposed.  The  predicted  forces  are  compared  to 
those  measured  in  existing  experiments. 
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2  UNIFORM  INFLOW 
2.1  Partial  wake  alignment 

This  model  was  developed  first  for  the  prediction  of  the  propeller  performance  in  steady  or  unsteady  flows 
by  using  a  VLM  approach  (Kerwin  k  Lee  1978).  The  corresponding  code  is  called  PUF-2.  The  trailing  wake 
is  divided  into  a  transition  and  an  ultimate  wake  region.  In  the  transition  region  the  trailing  vorticity  lines  are 
assumed  to  converge  to  either  the  tip  or  the  hub  roll-up  points  depending  on  the  sign  of  the  vorticity  strength1. 
The  roll-up  points  are  also  the  starting  points  for  the  concentrated  tip  or  hub  vortices  in  the  ultimate  wake 
region.  The  pitch  angle,  /?x ,  at  the  outer  extremity  of  the  transition  region  is  taken  by  the  following  relationship: 

a  =  +  hi 

2 

where  Pup  is  the  hydrodynamic  pitch  angle  of  the  undisturbed  inflow  and  <£0.7  is  the  geometric  pitch  angle  at 
r  =  0.7  R. 

The  pitch  angle,  Pw ,  of  the  concentrated  tip  vortex  in  the  ultimate  wake  region  is  taken  from  the  following 
relationship: 


Pw  =  1.15  /3t  (2) 

The  extent  of  the  transition  wake  is  determined  by  the  blade  angle  between  the  tip  of  the  blade  and  the 
beginning  of  the  ultimate  wake.  It  is  usually  taken  equal  to  90°.  The  radius  of  the  ultimate  wake  is  user  specified 
and  taken  equal  to  0.83  of  the  propeller  radius2. 

This  wake  model  was  also  applied  in  the  case  of  cavitating  flow  (Lee  1979).  The  corresponding  code  is  called 
PUF-3.  It  was  found  that  this  wake  model  was  resulting  to  cavity  sources  in  the  trailing  wakejm  the  case  of 
super-cavitating  flows)  the  total  value  of  which  was  inconsistent  with  the  value  of  the  cavity  volume  velocity 
(Kinnas  1983).  This  inconsistency  was  resolved  by  not  allowing  for  either  roll-up  or  contraction  of  the  trailing 
wake  lines  (Kinnas  1983).  The  revised  code  is  called  PUF-3A3.  The  corresponding  wake  model  is  shown  in  Figure 
1.  Note  that  the  transition  region  has  been  doubled  (from  a  blade  angle  of  90°  to  180°)  to  accommodate  longer 
super-cavities  in  one  of  the  most  recent  versions,  HPUF-3AL  (Pyo  k  Kinnas  1996). 

1  This  technique  is  accounting  for  the  increase  in  the  strength  of  the  tip  vortex  in  a  semi-empirical  manner  without  actually 
treating  the  mechanism  of  the  “wake  roll-up” 

2  This  value  was  measured  over  a  wide  range  of  propeller  geometries  and  advance  ratios  (Kerwin  &  Lee  1978). 

3  This  version  also  included  the  leading  edge  correction  (Kinnas  1991,  Kerwin  et  al  1986). 
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2.2  Wake  alignment  without  roll-up 

Equation  (1)  provides  partial  dependency  of  the  trailing  wake  on  the  incoming  flow.  An  iterative  procedure 
that  aligns  the  trailing  wake  geometry  with  the  resulting  flow  was  developed  in  (Greeley  &  Kerwin  1982).  The 
corresponding  code  is  called  PSF-2.  The  wake  is  again  divided  into  the  transition  and  the  ultimate  region.  The 
trailing  vortex  lines  are  aligned  with  the  actual  (computed)  flow  at  the  start  (blade  trailing  edge*)  and  at  the  end 
of  the  transition  wake  region,  while  interpolated  in  between.  Only  the  axial  and  tangential  velocities  induced 
by  the  propeller  are  being  used  in  this  procedure.  The  radial  location  of  the  vortex  lines  is  determined  in  terms 
of  the  contraction  of  the  transition  wake  at  the  tip  and  the  ultimate  wake  radius,  both  of  which  are  inputs 
taken  from  experimental  observations.  In  this  way  the  wake  roll-up  is  suppressed.  The  ultimate  wake  pitch  is 
determined  by  using  the  method  developed  in  (Loukakis  1971)  where  the  ultimate  wake  is  assumed  to  extent 
infinitely  upstream  and  downstream.  This  wake  model  is  depicted  in  Figure  2. 

This  method  was  also  applied  in  the  case  of  the  representation  of  the  blade  with  a  boundary  element  (panel) 
method  for  steady  (Lee  1987,  Kerwin  et  al  1987)  or  unsteady  flows  (Hsin  1990,  Kinnas  &  Hsin  1992).  The  names 
of  the  corresponding  codes  are  PSF-10  and  PUF-10,  respectively.  In  the  case  of  non-axisymmetric  inflows  the 
trailing  wake  shape  is  assumed  to  be  the  same  at  all  blade  angles  and  is  determined  from  alignment  with  the 
circumferentially  averaged  inflow. 

2.3  Improved  wake  model  in  PUF-3A 

The  model  described  in  Section  2.2  was  recently  implemented  in  PUF-3A  (Pyo  &  Kinnas  1996).  The 
corresponding  code  is  called  HPUF-3AL  and  also  includes  the  effects  of  the  hub  on  blade  cavitation  (Kinnas 
1996).  In  this  wake  model  the  ultimate  wake  radius  was  first  taken  equal  to  0.83  of  the  propeller  radius.  However, 
to  avoid  the  inconsistency  in  PUF-3A  (as  described  in  Section  2.1)  an  ultimate  wake  radius  equal-io  the  propeller 
radius  is  considered.  The  original  wake  (default)  and  the  aligned  wake  with  two  different  ultimate  wake  radii  are 
shown  in  Figure  3  for  propeller  DTMB  4661  (Boswell  et  al  1984). 

The  effect  of  the  wake  model  on  the  predicted  cavity  shape  can  be  seen  in  Figure  4.  Note  that  the  cavity 
extent  and  volume  are  not  affected  appreciably.  The  effect  of  the  wake  model  on  the  predicted  blade  forces  (in 
steady  wetted  flow)  is  shown  in  Table  1.  Note  that  the  current  model  underpredicts  the  forces  substantially 
when  the  old  wake  model  is  employed  in  HPUF-3AL.  In  addition,  the  forces  are  not  considerably  altered  when 
an  ultimate  wake  radius  of  1.0  instead  of  0.83  of  the  propeller  radius  is  considered. 
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WAKE  MODEL 

Fx 

Fy 

Fz 

WAKE  MODEL 

Qx 

Qy 

. Qz . 1 

Old  model 

-0.07569 

IAMI 

Old  model 

0.01651 

0.00020 

MiWiWiZU 

HI 

Aligned  (RULT=1) 

0.01728 

0.00031 

0.02719 

Aligned  (RULT=0.83) 

mmaem 

0.04905 

Aligned  (RULT=0.83) 

0.01725 

0.00036 

0.02720 

PSF-10 

|  -0.07852 

PSF-10 

0.01724 

fjgjMSlB; 

Table  1:  The  steady  fully  wetted  forces  per  blade.  Predicted  by  HPUF-3AL  with  default  and  aligned  wake 
models,  and  PSF-10  (with  aligned  wake  model).  Propeller  DTMB  4661,  Js  =  0.8. 


2.4  Full  wake  alignment  (with  roll-up) 

As  mentioned  in  Section  2.2,  the  iterative  lifting  surface  method  for  the  wake  alignment  “artificially”  sup¬ 
presses  the  tip  vortex  roll-up. 

The  contraction  of  the  wake  (which  is  due  to  the  mechanism  of  the  roll-up)  and  the  ultimate  wake  radius  are 
parameters  determined  from  experimental  observations.  The  full  wake  alignment  for  propellers  in  non-uniform 
inflows  was  first  treated  by  (Keenan  1989).  He  applied  a  vortex  lattice  method  (VLM)  on  both  blade  and  wake 
and  was  able  to  predict  the  gross  characteristics  of  the  wake  geometry.  However,  due  to  inherent  numerical 
instabilities  when  finer  discretization  is  used  or  during  late  stages  of  the  wake  evolution,  this  method  cannot 
predict  the  detailed  structure  of  the  wake  roll-up. 

In  (Pyo  1995)  and  (Pyo  Sc  Kinnas  1997),  an  iterative  process  for  the  complete  wake  alignment,  including 
roll-up,  was  introduced  and  implemented  in  the  propeller  panel  method.  This  method  was  also  coupled  with  the 
flow  adapted  grid,  FLAG,  introduced  in  (Kinnas  et  al  1993).  The  corresponding  code  is  called  PSF-10FLAG. 
The  location  where  the  tip  vortex  detaches  from  the  blade  is  defined  by  the  ratio  AS/ AT,  as  depicted  in  Figure 
5.  Panels  with  bi-quadratic  dipole  strength  are  utilized  in  order  to  model  the  wake  surface  and  in  particular 
the  wake  roll-up  region.  The  method  is  capable  of  predicting  the  detailed  structure  of  the  wake  roll-up  in  three 
dimensions  as  shown  in  Figure  8.  The  contraction  of  the  wake  and  the  ultimate  wake  radius  are* now  evaluated 
by  the  method. 

It  should  be  noted  that  for  some  blade  geometries  and  inflow  conditions  the  predicted  wake  sheet  roll-up 
may  be  sensitive  to  the  number  of  panels  in  the  wake,  despite  the  high-order  dipole  panel  representation.  This 
has  been  attributed  to  the  inadequacy  of  coarser  grids  to  capture  the  slope  of  the  blade  loading  at  the  tip,  which 
in  turn  has  been  found  to  be  very  critical  in  evaluating  the  initial  stage  of  the  roll-up.  In  order  to  get  sufficient 
representation  of  the  roll-up  during  its  initial  stage  a  minimum  number  of  panels  has  been  found  to  be  required  in 
the  wake  (usually,  40  spanwise  and  30  streamwise  panels).  In  addition,  the  wake  shape  and  induced  velocities  are 
re-splined  after  each  roll-up  iteration.  These  changes  (described  in  more  detail  in  Appendix  A)  were  incorporated 
in  PSF-10FLAG  and  were  found  to  render  reliable  results  over  a  wide  range  of  blade  geometries  and  advance 
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ratios. 


The  location  of  the  tip  detachment  point  can  affect  the  results  at  the  tip  substantially.  This  is  shown  for 
example  in  Figure  6  where  the  pressure  distribution  on  the  blades  of  DTMB  5168  is  shown  for  three  different 
locations  of  the  tip  vortex  detachment.  Note  that  the  pressure  distribution  everywhere  on  the  blade  (except 
at  the  tip)  is  unaffected'  by  the  detachment  point.  In  this  particular  case  AS  =  1.20 AT  appears  to  be  the 
detachment  point  with  the  smallest  region  of  singular  regions  at  the  tip.  The  corresponding  blade  and  rolled-up 
wake  grids  are  shown  in  Figure  7.  More  recently,  a  wake  roll-up  process  based  on  the  de-singularization  of  the  tip 
vortex  (Krasny  1987)  was  developed  by  (Ramsey  1996).  In  the  case  of  wings  a  study  of  the  tip  vortex  detachment 
location  on  the  tip  flow  is  described  in  Appendix  B.  Finally,  the  prediction  of  the  viscous  flow  inside  the  core  of 
a  tip  vortex  is  addressed  in  Appendix  C. 

2.5  Comparison  with  experiments 

The  predicted  wake  shape  by  PSF-10FLAG  for  DTMB  4119  is  shown  in  Figure  8  to  be  very  close  to  measured 
(Jessup  1989). 

The  predicted  open  water  characteristics  by  different  methods  are  shown  together  with  the  measured  in 
Figure  9.  The  geometry  and  the  open  water  measurements  for  DTMB  5168  were  provided  by  (Jessup  1996). 
The  geometry  is  also  given  in  Appendix  D.  A  friction  coefficient  of  0.0035  has  been  used  on  all  methods.  All 
methods  seem  to  give  very  good  predictions  around  design  0.8  <  Js  <  1.2.  HPUF-3A-INC  (same  as  HPUF-3AL 
for  zero  shaft  inclination)  appears  to  underpredict  Kt  for  Js  <  0.8  and  Kq  for  Js  >  1.2.  PSF-104  appears  to 
underpredict  the  torque  for  Js  >  1.2.  Finally,  PSF-10FLAG,  appears  to  predict  forces  closer  to  experiment  for 
a  wide  range  of  0.6  <  Js  <  1.6.  The  improvement  in  the  Kq  predicted  by  PSF-10FLAG  may  be  explained  by 
the  fact  that  the  detachment  point  with  the  “smallest  singular  pressure  region”  was  used  in  the  calculation,  and 
that  the  moment  of  the  tip  pressures  with  respect  to  the  propeller  axis  can  be  appreciable.  More  investigation, 
including  systematic  observations  of  the  location  of  the  tip  vortex  detachment  point,  is  needed  to  confirm  our 
conjecture  regarding  the  improved  PSF-10FLAG  predictions  at  high  J’s. 

The  rolled-up  wake  geometry  predicted  by  PSF-10FLAG  for  Js  =  0.6  is  shown  in  Figure  10.  It  should  be 
noted  that  the  wake  is  assumed  to  start  at  the  propeller’s  tip  (AS /AT  =1)  without  allowing  for  leading  edge 
vortex  separation  which  may  be  present  at  this  low  advance  ratio.  In  regards  to  the  stability  of  our  method, 
as  also  mentioned  in  Section  2.4,  a  sufficient  number  of  panels  is  needed  in  the  spanwise  direction  in  order  to 
produce  stable  roll-ups. 

4This  version  of  PSF-10  uses  a  wake  alignment  scheme  similar  to  that  described  in  Section  2.2. 
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3  INCLINED  INFLOW 
3.1  Trailing  wake  geometry 


A  propeller  subject  to  uniform  inflow  with  inclination  angle  a  is  shown  in  Figure  11.  This  is  equivalent5 
to  having  the  propeller  subject  to  uniform  inflow  with  the  shaft  being  inclined  at  an  angle  a  (positive  when 
the  downstream  part  of  the  shaft  is  at  a  larger  depth  than  the  upstream  part).  The  total  inflow,  Vin  may  be 
decomposed  into  an  axial  component,  Va  =  V/jv  cos  a,  and  a  vertical  component,  Vy  =  Vim  sin  a,  as  also  shown 
in  Figure  11.  Subsequently,  Vy  can  be  decomposed  into  a  radial  component,  Vr,  and  a  tangential  component, 
Vt,  given  as: 


Vr  (6)  =  Vy  cos  6  =  Va  tan  a  cos  9  (3) 

Vt{9)  =  —  Vy  sin#  =  —Va  tanasin#  (4) 

The  unsteady  propeller  flow  in  the  presence  of  the  uniform  axial  and  the  once  per  revolution  varying  radial 
and  tangential  inflows  must  be  solved.  As  already  mentioned  in  the  introduction  the  alignment  of  the  trailing 
wake  for  each  blade  angle  is  required  in  order  to  predict  the  blade  forces  accurately.  Ideally,  the  trailing  wake 
geometry  could  be  determined  at  each  blade  angle  by  alignment  with  the  local  (time  dependent)  flow.  Since  this 
process  is  expected  to  be  prohibitively  computer  intensive,  a  more  efficient  procedure  is  proposed  in  this  paper. 
This  procedure  is  depicted  in  Figures  12  and  13.  The  trailing  wake  geometry  is  determined  as  the  superposition 
of  (a)  the  trailing  wake  geometry  for  uniform  axial  inflow  with  zero  inclination,  and  (b)  an  angular  shift  of  the 
wake  geometry  equal  to  the  angle  of  inclination.  As  a  result  of  this  superposition  the  trailing  wake  geometry  will 
change  with  blade  angle  as  shown  is  Figure  13.  This  approximation  of  the  wake  geometry  is  expected  to  be  as 
accurate  in  inclined  flows  as  the  approximations  involved  in  the  method  of  (Greeley  &  Kerwin  1932)  in  straight 
inflows. 

In  Figure  11  the  total  velocity  vector  VQ  (average  of  the  flow  velocities  on  the  two  sides  of  the  wake  sheet)  on 
the  trailing  wake  surface  is  shown  with  respect  to  a  coordinates  system  rotating  with  the  propeller.  The  axial, 
radial,  and  tangential  components  of  the  total  velocity  (shown  in  Figure  11)  are  given  as  follows: 


vA 
y  0 

-  VA  +  ua 

(5) 

vR 

V  0 

=  vR{0) 

(6) 

VT 

V  0 

=  Vr(0)  +  ut+u>r 

(7) 

5  We  ignore  the  wake  defect  in  the  shadow  of  the  inclined  shaft. 
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where  ua  and  ut  are  the  axial  and  the  circumferential  components  of  the  propeller  induced  perturbation  velocity, 
respectively.  These  velocities  are  evaluated  with  the  propeller  subject  only  to  the  uniform  axial  inflow  (i.e.  at 
zero  inclination  angle).  The  method  described  in  Section  2.2  is  used  at  this  stage.  Note  that  the  induced  radial 
perturbation  velocity  is  not  included  in  V0,  since  this  component  is  not  utilized  in  the  wake  alignment  process, 
as  explained  in  an  earlier  section. 

The  geometry  of  the  transition  wake  is  determined  by  starting  at  the  trailing  edge  at  each  radial  position  and 
then  advancing  in  constant  increments  AO  in  the  circumferential  direction.  Usually  A 6  =  6°.  The  increments  in 
terms  of  the  edges  of  the  panels  in  the  trailing  wake  are  given  from  the  following  equations.  This  is  also  depicted 
in  Figure  11. 


Ax  =  VcaA t  =  VAysUa  ( Js 

A0\ 

180  ) 

(8) 

Ay  =  ARcosO  —  ATsinO 

(9) 

A z  =  Aiisin  6  +  AT  cos  6 

(10) 

A6  _  Js  A 6 
u  ~  Vs  180 

(11) 

ar  =  ^ !  (*®) 

(12) 

at  =  nr* =[*+•'+;,] 

("&)  ' 

(13) 

where  a  is  the  inclination  angle,  0  is  the  angular  position  in  the  wake  of  the  blade,  Js  is  the  advance  coefficient, 
w  is  a  propeller  angular  velocity.  In  the  ultimate  wake  region,  the  pitch  is  calculated  in  the  same  manner  as  in 
the  uniform  inflow  case.  The  aligned  wake  geometries  at  three  angular  positions  are  also  shown  in  Figure  14. 

This  wake  model  has  been  implemented  in  HPUF-3AL  and  the  new  version  is  called  HPUF-3A-INC.  The 
influence  coefficients  of  the  trailing  wake  vortex  horse-shoes  on  the  key  blade  are  evaluated  at  each  blade  angle 
(60  in  total).  As  a  result  the  computing  time  is  increasing  considerably  (about  five  times)  when  the  inclined 
wake  geometry  is  included.  A  20  x  18  HPUF-3A-INC  run  on  a  DEC  Alpha  600-5/266  takes  about  15  minutes  of 
CPU  time  when  the  inclined  wake  geometry  is  included. 

3.2  Results 

Results  are  presented  next  for  propeller  DTMB  4661  (Boswell  et  al  1984)  the  geometry  of  which  is  given  in 
Appendix  D. 
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First  the  effect  of  the  inclined  trailing  wake  geometry  on  the  blade  forces  is  shown  in  Figure  15.  Note  the 
considerable  difference  when  the  inclined  wake  geometry  is  included  in  the  calculation.  A  similar  effect  was  also 
reported  in  (Kerwin  1979). 

Figures  16  and  17  show  the  effect  of  the  inclined  wake  geometry  on  the  predicted  cavity  shapes.  The  inclined 
wake  model  predicts  a  larger  cavity  extent. 

Figures  18  and  19  show  the  predicted  cavity  volume  for:  (a)  straight  trailing  wake  and  by  considering  only 
the  axial  inflow  component,  (b)  same  as  a)  by  also  adding  the  radial  component,  (c)  same  as  b)  by  also  adding 
the  tangential  component  and  (d)  inclined  trailing  wake  with  all  inflow  components.  Note  the  expected  negligible 
effect  of  the  radial  component  on  the  predicted  cavity  volume.  The  inclusion  of  the  tangential  component  in  the 
inflow  changes  the  value  and  the  angular  location  of  the  maximum  cavity  volume  drastically.  This  is  due  to  an 
increase  of  the  local  angle  of  attack  at  90°  and  a  decrease  of  the  local  flow  angle  at  270°,  when  the  tangential 
inflow  component  is  considered.  Finally,  the  cavity  volume  increases  even  more  when  the  inclined  wake  model  is 
included.  This  increase  is  consistent  with  the  already  shown  increase  in  the  blade  forces. 

Figures  20  and  21  show  the  amplitudes  of  the  first  harmonic  of  the  forces  acting  on  the  blades  of  propeller 
DTMB  4661  for  two  inclination  angles  and  over  a  range  of  advance  ratios.  The  measured  forces  (Boswell  et  al 
1984)  are  shown  together  with  those  predicted  by  the  present  method  and  the  method  in  (Kerwin  1979)  where 
the  inclined  wake  geometry  is  considered  at  blade  angles  of  0°,  90°,  180°  and  270°.  Note  the  fair  agreement  of 
the  predictions  from  both  methods  with  the  measurements.  The  present  method,  however,  incorporates  a  more 
complete  inclined  wake  model  (which  may  be  necessary  for  propeller  geometries  or  inclinations  other  than  those 
presented  in  Figures  20  and  21),  and  most  importantly,  includes  the  modeling  of  cavitation. 

4  CONCLUSIONS  - 

Existing  trailing  wake  models  are  reviewed  and  the  predicted  open  water  characteristics  are  compared  again 
each  other  and  against  experimental  measurements.  While  the  more  sophisticated  models  can  predict  wake 
shapes  and  pressure  distributions  at  the  blade  tip  more  accurately  than  the  simpler,  all  models  seem  to  predict 
the  forces  within  acceptable  accuracy,  especially  in  the  range  of  design  J. 

One  of  the  simpler  models  is  extended  in  the  case  of  cavitating  flow  in  inclined  inflow.  The  geometry  of  the 
inclined  wake  is  calculated  at  all  blade  angles.  The  effect  of  the  inclined  wake  on  the  predicted  cavity  volume 
and  unsteady  blade  forces  is  found  to  be  significant.  The  predicted  forces  compare  well  to  those  measured  in 
inclined  shaft  propeller  experiments. 
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Figure  1:  Partially  aligned  wake  model  in  PUF-3A  (original  version) 


Figure  2:  Aligned  wake  model  in  PSF-2. 


15 


m 


Figure  7:  Fully  aligned  wake  model  in  PSF-10FLAG.  Propeller  DTMB  5168  at  Js  =  1.0. 
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Figure  8:  Radial  location  of  the  tip  vortex  for  propeller  DTMB  4119.  From  Pyo  and  Kinnas  (1997). 


Figure  9:  Kt$z10Kq  vs.  J s  for  propeller  DTMB  5168;  measured  and  predicted  by  different  methods. 
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Figure  10:  Wake  geometry  aligned  using  PSF-10FLAG.  Propeller  DTMB  5168  at  Js  =  0.6. 


Figure  11:  A  propeller  subject  to  uniform  inflow  with  inclination  a  and  the  geometry  of  the  inclined  trailing 
wake  geometry. 
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Figure  14:  The  trailing  wake  geometry  aligned  with  the  inclined  flow  shown  at  three  angular  positions; 
a  =  20°. 


Figure  15:  The  first  harmonic  of  the  forces  and  moments  acting  on  one  blade  for  10°  shaft  inclination; 
predicted  by  the  present  method  with  straight  and  inclined  trailing  wake.  Propeller  DTMB  4661,  fully 
wetted  flow. 
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Figure  16:  The  cavity  shape  at  four  angular  positions  for  10°  shaft  inclination;  predicted  by  HPUF-3AL 
(straight  aligned  wake)  and  HPUF-3A-INC  (inclined  aligned  wake).  Propeller  DTMB  4661,  Js  =  1.0,  a  = 
0.2,  F„  =  1.165. 


Figure  17:  The  cavity  shape  at  four  angular  positions  for  20°  shaft  inclination;  predicted  by  HPUF-3AL 
(straight  aligned  wake)  and  HPUF-3A-INC  (inclined  aligned  wake).  Propeller  DTMB  4661,  Js  =  1.0,  a  = 


0.2,  Fn  =  1.165. 
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Figure  18:  The  cavity  volume  as  a  function  of  the- blade  angle  for  10°  shaft  inclination;  predicted  by  HPUF- 
3A-INC  for  different  inflows  and  trailing  wake  geometries.  Propeller  DTMB  4661,  Js  =  1.0,  cr  =  0.2,  Fn  = 
1.165. 
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Figure  19:  The  cavity  volume  as  a  function  of  the  blade  angle  for  20°  shaft  inclination;  predicted  by  HPUF- 
3A-INC  for  different  inflows  and  trailing  wake  geometries.  Propeller  DTMB  4661,  Js  —  1.0,  a  =  0.2 ,F„  = 
1.165. 
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Figure  20:  The  first  harmonic  of  the  forces  and  moments  acting  on  one  blade  for  10°  shaft  inclination; 
measured  and  predicted  by  the  present  (HPUF-3A-INC)  and  a  previous  method  (PUF-2IS).  Propeller  DTMB 
4661  in  fully  wetted  flow. 


Figure  21:  The  first  harmonic  of  the  forces  and  moments  acting  on  one  blade  for  20°  shaft  inclination; 
predicted  by  the  present  (HPUF-3A-INC)  and  a  previous  method  (PUF-2IS).  Propeller  DTMB  4661  in  fully 
wetted  flow. 
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A  The  algorithm  of  the  wake  alignment  model 

III.  ALIGNED  WAKE  MODEL  WITH  ROLL-UP 


•  PROCEDURE 

(I)  FIRST  ITERATION 


FLAG  wl  0  angle  Straight  wake 
(2)  SECOND  &  MORE  ITERATIONS 
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(b)  IMPROVEMENT  ON  LOADING  &  FORCE  CALCULATION 


-}-  DIFFERENCE  BETWEEN  BEFORE  I.P.C.  AND 
AFTER  I.P.C.  BECOMES  SMALLER. 

-}•  FAST  CONVERGENCE  ON  I.P.C. 
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IV.  NUMERICAL  RESULTS 


•  CIRCULAR  PLANFORM  WING  . 

(a)  MORE  STABLE  &  RELIABLE  PREDICTION  ON  THE 
SHAPE  OF  TRAILING  WAKE  SHEET  ROLL-UP. 
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ON  THE  BLADE 

(a)  B-SPLINE  CURVE  WITH  .  GIVEN  ANGLE  AT  THE 

TRAILING  EDGE  &  NORMAL  TO  L.E..  ... 

(b)  HALF  COSINE  SPACING  ALONG  SPANWISE 
DIRECTION  &  COSINE  SPACING  ALONG  CHORDWISE 
DIRECTION. 

(c)  TWO  SAME  SIZE  PANELS  AT  THE  TAILING  EDGF 

(d)  GET  TRAILING  EDGE  ANGLE  FROM  VELOCITY  OF  THE 
T.E.  PANEL  ON  THE  BLADE. 

IN  THE  WAKE 

(a)  HIGH  ORDER  PANEL  METHOD  WITH  REDISCRETIZATION 

(b)  FIXED  NUMBER  OF  PANELS  (20  X  30) 

(c)  HALF  COSINE  SPACING  ALONG  STREAMWISE 
DIRECTION  &  COSINE  SPACING  ALONG  SPANWISE 
DIRECTION. 


B  Effect  of  detachment  location  on  tip  flow  for  wings 

DEFINITION  OF  TIP  VORTEX  DETACHMENT 

POINTS 


S  at  x  :  AS  =~r  AT 


1)  Effect  of  detachment  point  on  wake  shape 
at  the  tip/TE  for  thin  wings,  (t/c  =  2%,  6°  angle 
of  attack) 

.m 

•  l.oo  - 
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1)  Effect  of  detachment  point  on  wake  shape 
at  the  tip/TE  for  thick  wings,  (t/c  =  20%,  6° 

angle  of  attack)  ’ 

-•  * 

1.00 
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1)  Effect  of  detachment  point 'on  wake  shape 
at  the  tip/TE  for  thick  wings,  (t/c  =  20%,  10° 
angle  of  attack)  ’ 


2)  Effect  of  detachment  point  on  pressure 
distribution  (t/c  =  20%,  6°  angle  of  attack)  , 


PRESSURE  SIDE 


SUCTION  SIDE 


PRESSURE  SIDE 


:  ‘-V 


C  Local  viscous  flow  inside  of  a  tip  vortex 

V.  VISCOUS  EFFECT  INSIDE  THE  CORE 
1.  FORMULATION 

•  ASSUMPTION  :  SLENDER  CORE  &  STEADY  &  AXIS-YMMETRY 

•  GOVERNING  EQUATION 


d(rv)  d(rv) 
+  w-  v  J 


dx 


dr 


du  du 
u— +  w— 
ox  dr 


dw  dw 
u— — h  w 


v 


dx  dr  r 
du  dw  w 

dx+  dr  +  r 


•  USE  FINITE  DIFFERENCE  METHOD  (PARABOLIC  EQUATION) 
CENTRAL  DIFFERENCE  ALONG  r  DIRECTION  ‘ 

BACKWARD  DIFFERENCE  ALONG  x  DIRECTION 


INITIAL  CONDITION 


P  ( x ,  r ) :  PRESSURE  OUTSIDE  CORE  (  INVISCID  FLOW) 
L :  CORE  DIAMETER 
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2.  EXAMPLE  (A  TRAILING  VORTEX) 

♦  CONSTANT  CORE  RADIUS  ALONG  THE  CORE  CENTER  LINE. 

m  m 

♦  CONSTANT  PRESSURE  OUTSIDE  THE  CORE. 

•  INITIAL  CONDITION  : 

x0  =  0.25,  R(x0)  =  6 

u0  =  1  —  0.25e“r  ,  v0  =  (0.5  —  0.5 e~r2)/r 

•  AT  x-x0  :  CORE  RADIUS  (R0) 

VELOCITIES  INSIDE  THE  CORE  ( u0  ,  lb  ) 

PRESSURE  OUTSIDE  THE  CORE(P(x)) 
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CIRCUMFERENTIAL  VELOCITY  IN  A  TRAILING  VORTEX 


AXIAL  VELOCITY  IN  A  TRAILING  VORTEX 
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D  Propeller  geometries 


Number  of  Blades 
Hub/Diameter  Ratio 
Section  Meanline 
Section  Thickness  Form 
Design  J5 


5 

0.3 

NACA  a=0.8 
NACA66  (Modified) 
1.038 


pza 

skew0 

l!£21 

■wra-i 

0.0091 

2.985 

0.178 

0.0000 

[Iglggglj 

1.296 

0.0103 

0.210 

lililiHil 

■BBKBiij 

1.480 

0.0103 

4.810 

liEOT 

0.0290 

eei 

■if.W 

0.0267 

0.0226 

1.566 

0.0103 

1.498 

0.0103 

11.895 

0.406 

0.0209 

■iKOEftM 

1.381 

0.0103 

15.410 

0.409 

0.0151 

BTIiYKRl 

WmUM 

0.0102 

0.387 

0.0122 

i«EESl 

0.0103 

0.326 

0.0094 

0.0091 

21.876 

0.000 

0.0000 

0.0000 

Table  2:  The  geometry  of  the  propeller  DTMB  4661. 
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Number  of  Blades 
Hub/Diameter  Ratio 
Section  Meanline 
Section  Thickness  Form 
Design  Js 
Cavitation  Number 
Froude  Number 


5 

0.3 

NACA  a=0.8 
NACA66  (Modified) 
1.0 
2.0 
8.86 


r/R 

P/D 

%m  /  D 

skew0 

c/D 

f/c 

t/D 

0.30 

1.087 

0.0111 

6.280 

0.200 

-0.0380 

0.0539 

0.35 

1.245 

-0.0069 

-0.754 

0.237 

-0.0086 

0.0422 

0.40 

1.365 

-0.0237 

-4.824 

0.276 

0.0090 

0.0315 

0.45 

1.458 

-0.0358 

-7.336 

0.308 

0.0178 

0.0274 

0.50 

1.541 

-0.0437 

-8.865 

0.334 

0.0279 

0.0252 

0.60 

1.673 

-0.0483 

-9.838 

0.385 

0.0366 

0.0226 

0.70 

1.633 

-0.0419 

-8.108 

0.435 

0.0332 

0.0021 

0.80 

1.502 

-0.0302 

-3.784 

0.475 

0.0247 

0.0019 

0.85 

1.428 

-0.0234 

-0.493 

0.482 

0.0186 

0.0186 

0.90 

1.335 

-0.0165 

3.784 

0.465 

0.0119 

0.0203 

0.93 

1.259 

-0.0123 

6.922 

0.431 

0.0079 

0.0227 

0.95 

1.189 

-0.0096 

9.297 

0.390 

0.0054 

0.0229 

0.97 

1.097 

-0.0068 

11.927 

0.325 

0.0030 

0.0205 

0.99 

0.974 

-0.0039 

14.834 

0.206 

0.0009 

0.0131 

1.00 

0.900 

-0.0025 

16.400 

0.000 

0.0000 

0.0000 

Table  3:  The  geometry  of  the  propeller  DTMB  5168. 


